Effects of local oxygen distortions on electronic structures of Na^CoC^ 
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By using pseudopotential method with local spin density functional approximation, the electronic 
band structures of Na^CoC^ are calculated for x — 0.25, 0.5, 0.75, and x = 1 in the presence of 
the structure relaxations. As increasing Na content, the hybridization between cobalt and oxygen 
orbitals is decreased, and a phase transition is predicted from a wide-band ferromagnetic to a 
narrow band paramagnetic metals. The itinerant ferromagnetism is strongly suppressed by the local 
distortions of the oxygens around the cobalts. Moreover, straining the C0O2 layers corresponding 
to the hydrated superconductor Nao.35CoO2T.3H2O strongly enhances both the hybridization and 
ferromagnetism. 

PACS numbers: 71.20-b, 74.25. Jb, 74.70.Ad 
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Since the discovery of high critical temperature su- 
perconductivity in layered copper oxide compounds, the 
search for layered superconductors without copper atoms 
has been regarded as an important route to understand 
the unknown mechanism behind the highest supercon- 
ducting transition temperatures. Recently, Takada and 
coworkers pj have discovered a novel superconductor 
Nao.35CoO2T.3H2O with T c = 5K, which is synthesized 
when a sodium cobalt oxide Nao.7sCo0 2 is hydrated to 
reduce its sodium content and intercalated with water 
molecules to further separate the C0O2 layers. Based on 
the analogy of its layered transition metal oxide struc- 
ture, it has been speculated that the displayed super- 
conductivity may be closely related to that of the high 
T c cuprates iQlQ ■ Along this line, several theoretical 
analysis 0, U Q have been carried out in terms of res- 
onant valence bond states on a strongly correlated two- 
dimensional t-J model on triangular lattices. 

However, unlike the parent compounds of high T c 
cuprates as an antiferromagnetic Mott insulator, the par- 
ent compounds Na x Co0 2 (x < 1) have been known 
for their anomalous large thermoelectric power and low 
electronic resistivity Na^Co02 consists of a two- 

dimensional triangular lattice of cobalt ions formed by 
a network of edge-sharing C0O6 octahedra, separated by 
layers of sodium ions. Band structure calculations us- 
ing a general potential linearized augmented plane wave 
method 7] for Na .5CoO 2 indicates that it is a ferromag- 
netic half metal with little hybridization between the oxy- 
gen and cobalt orbitals, while no distinct magnetic order- 
ing has been observed experimentally 0, Q . But for the 
samples of Nao.7sCo02, a magnetic phase transition as- 
sociated with the antiferromagnetic interactions has been 
confirmed below T = 22K, being accompanied by a spe- 
cific heat jump, a small magnetization, and a kink in the 
temperature dependence of the resistivity |lOl Hlj . 

Moreover, it has been reported that the crystal 
structure of the superconducting Nao.35Co0 2 T.3H 2 
compound only differs from the parent compound 
Nao.75Co02 by the intercalation of water molecules to en- 
large the Co0 2 layers and to reduce the sodium charges 



P, . It is thus of great interest to scrutinize the elec- 
tronic band structures of the parent compound Na K Co02 
for different sodium contents and to study the electronic 
structure with a large C0O2 layers as the same as that 
of the hydrated compound. 

In this paper, by using pseudopotential method with 
local spin density functional approximation 0, 0, ^| , 
we present the first principles electronic structure calcu- 
lations on Na K Co0 2 for different Na content. The elec- 
tronic structure is found to be highly two dimensional, 
giving rise to a good metal for x = 0.25, 0.5, 0.75 and a 
semiconductor for x = 1. In particular, at x = 0.25 it is 
an itinerant ferromagnetic metal with large hybridization 
between cobalt and oxygen orbitals. As increasing the Na 
content, both the ferromagnetism and hybridization are 
decreasing. When x = 0.5, only tiny ferromagnetism is 
left, and the hybridization vanishes. A phase transition 
is thus yielded from the wide band ferromagnetic metal 
to a narrow band paramagnetic metal by changing the 
Na content. The ferromagnetism is totally absent and a 
narrow conduction band appears below the Fermi energy 
when x = 0.75. We find that such a strong suppression 
of the ferromagnetism is mainly originated from the local 
distortions of the oxygen ions around the cobalt atoms. 
However, the present treatments for the narrow band 
paramagnetic state may underestimate the tendency of 
the material towards local moment formation and possi- 
ble magnetic ordering. In addition, straining the sep- 
aration of the C0O2 layers hypothetically correspond- 
ing to the hydrated superconductor Nao.35CoO2T.3H2O, 
we find that both the hybridization and ferromagnetism 
are enhanced, and the corresponding superconductivity 
might be associated with a triplet (p-wave) pairing state 
in the presence of strongly ferromagnetic fluctuations 
[l5l IToil . analogous to those in UGe2, URhGe, and ZrZn 2 

The crystal structure used in the present numerical 
calculations is based on the neutron diffraction data re- 
ported by Balsys and Davis j^jj. Neutron data for 
Nao.74Co02 shows that it belongs to the hexagonal 
P63/mmc space group with lattice parameters a = 
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2.840A and c = 10.81lA. There are two Co0 2 sheets 
per unit cell. The triangles formed by the oxygen ions in 
the oxygen layers have two directions. We use A and B to 
represent these two oxygen layers. The oxygen packing in 
the structure is A ABB. The experimental results indicate 
that the oxygen atoms are at the sites (1/3, 2/3, ±0.088). 
The Na ions are intercalated with a trigonal prismatic site 
between the octahedral sheets of cobalt oxide. There are 
two available sites for the sodium ions: (2/3,1/3,1/4) 
and (0,0, 1/4). The sodium ions occupy these two sites 
unequally due to the influence of the Co atoms above 
and below in the second layer. The probability of find- 
ing sodium ions at the site (2/3, 1/3, 1/4) being approx- 
imately twice that on the other site (0, 0, 1/4) 20]. 

In order to take into account the different stoichiome- 
try x = 0.25, 0.5, and 0.75, an enlarged unit cell (2x2x1) 
of the crystal structure is introduced. The electronic 
Kohm-Sham orbitals for the momentum k-points in the 
Brillouin zone are expanded in plane waves with a cut- 
off of 24 Ryd. The PWSCF code is used in the 
calculations with the ultrasoft pseudopotential |22l and 
Perdew-Wang exchange correlation potential [14|, Us- 
ing this method, we can properly treat the atomic re- 
laxations and consider their relations with the electronic 
structures. We further assume the sodium ion setting on 
the site (2/3, 1/3, 1/4). In fact, in our framework the po- 
sition of sodium ions have little effect on the electronic 
structures, because the role played by sodium ions mainly 
donates charges to the C0O2 layers. More importantly, 
we have performed the structural relaxation to determine 
the optimized atomic positions. In particular, it is found 
that there are additional local distortions of the oxygen 
ions. While the earlier first principle calculations |7L Hq 
did not consider such local distortions. 

During the calculations, we have noticed that the lo- 
cal distortions of the oxygen ions strongly depend on 
the sodium content. As increasing the Na content, the 
oxygen distortions are large, and the itinerant ferromag- 
netism of Na^CoC^ is reduced significantly. To describe 
these local distortions, we define that Zi denotes the 
heights of the oxygen ions of each layer in their z-direction 
to the Co layer in the unit of lattice spacing c. All the 
oxygen layers (A and B) have the same magnitudes of 
the local distortions. The relaxed positions of the oxy- 
gen ions have been evaluated with respect to the sodium 
content and given in Table 1. We would like to point out 
that the case with sodium ions on the site (0,0,0) has 
similar relaxations for the oxygen ions. 

Table I. The oxygen positions and spin magnetizations 
per cobalt atom for the unrelaxed (m u ) and relaxed struc- 
tures (m r ) with respect to the sodium content (x). 



X 


z,(c) (1=1,2,3,4) 


m u (fi B ) 




x=0.25 


zi=0.0750; z 2 ,z 3 ,Z4=0.0787 


0.775 


0.340 


x=0.5 


zi,z 2 =0.0826; z 3 ,z 4 =0.0773 


0.534 


0.018 


x=0.75 


zi=0.0861; z 2 ,z 3 ,z 4 =0.0824 


0.263 





x=l 


z^0.0921 









The electronic band structures are displayed in Fig.l 
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FIG. 1: The electronic band structure for Nao.2sCo02. The 
Fermi energy Ef has been set at zero and denoted by the 
dotted lines. The left column corresponds to the spin-up elec- 
trons and the right column to the spin-down electrons. 



and Fig. 2 for x = 0.25 and x = 0.75, respectively. The 
energy scale is relative to the Fermi energy, and the den- 
sity of states is in atomic unit (Ryd) -1 . In Fig.l, the 
left column corresponds to the spin minority electronic 
structure, and the right column to the spin majority one. 
Both the spin minority and majority electronic structures 
are the same for Nao.75Co0 2 . The spin density of states 
(DOS) is calculated and delineated in Fig. 3 for the Na 
content x = 0.25, 0.5, 0.75, and 1. It is seen that a fer- 
romagnetic metal with the wide hybridized bands is ob- 
tained for x — 0.25. For x = 0.5, there is a very weak fer- 
romagnetism and little hybridization between the cobalt 
and oxygen orbitals. Moreover, a paramagnetic metal 
with a narrow band below the Fermi energy is found for 
x = 0.75. By adding electron charges to the Co0 2 lay- 
ers, the Fermi energy approaches to the top edge of the 
conduction band, and the difference between the spin mi- 
nority and spin majority is decreasing quickly. For the 
case of a; = 1, NaCo0 2 becomes a semiconductor, con- 
sistent with the earlier experimental result. However, in 
the earlier density functional calculations without the lo- 
cal distortions of oxygen ions a ferromagnetic half 
metallic state is obtained for all the sodium contents in 
x = 0.3, 0.5, and 0.7. In fact, we have also reproduced 
these results within the present calculation scheme. 

Moreover, the local spin density functional calculations 
can also self-consistently determine the spin magnetiza- 
tion of the ferromagnetically ordering state for low Na 
content (x < 0.5). In Table 1, the corresponding spin 
magnetizations per Co atom are given for both unrelaxed 
and relaxed structures. For x — 0.5, our result with the 
unrelaxed oxygen position is 0.534/Xg per cobalt atom. 
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FIG. 2: The electronic band structure for Nao.7sCo02. There 
are no differences between the spin-up and spin-down elec- 
trons. The Fermi energy Ep has been set at zero and denoted 
by the dotted lines. 



However, in the relaxed structure of material, it is the 
local oxygen distortions that reduce the spin magnetiza- 
tion to 0.018/^ B per cobalt ion. This result is consistent 
with the experimental observation |(| and the nominal 
valence of +3.5 of each cobalt ion in Nao.sCoC^. There- 
fore, it is expected that the material Nao.5Co02 is in 
the vicinity of the critical sodium content separating the 
weakly coupling wide band ferromagnetic metal from the 
strongly coupling narrow band paramagnetic metal. Par- 
ticularly, the Coulomb interaction in Na .75CoO2 may be 
larger than the conduction electron bandwidth, and the 
present spin density functional approximation may not 
be adequate to considering the possible antiferromagnetic 
instability. 

In order to exhibit the anisotropy of the charge 
distribution, we also calculate the local density of 
states (LDOS) at Fermi energy defined by n(Ep,r) = 
J2i l^'i( r )| 2 ^(^'F — which can provide some insights 
into the electron bonding around the individual atoms. 
The LDOS at the Fermi energy contour maps for the 
(001) surface are given in Fig. 4. For both Nao.7sCo02 
and Nao.25Co02, Co ions have a spherically symmetric 
density distribution around each site on the Co layer 
(Fig. 4b and 4e), while their local density distributions 
change from a down-triangular symmetry below Co layer 
(Fig. 4a and 4d) to an up-triangular symmetry above 
Co layer (Fig. 4c and 4f). The down- and up-triangles in 
the lower sodium content have more sharper edges than 
the higher sodium content, and the oxygen ions with a 
spherically symmetric distribution have higher density in 
low sodium content as well. There is a little charge den- 
sity around Na ions, reflecting that the role of Na atoms 
is mainly to donate charges to the C0O2 layer. These 
results clearly exhibit that a mixed metallic and covalent 
bonding are associated with the hybridization between 
the cobalt and oxygen orbitals, and that the compounds 
are clearly quasi-two dimensional systems. 

Finally, to understand the hydration effects of the 
compounds Na x CoC>2 (x < 1) on the electronic band 
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FIG. 3: The density of states of Naa;Co02 for the different 
sodium content with relaxed structure. The density of states 
is in atomic unit (Ryd) -1 . The up-spin DOS is shown above 
the horizontal line and the down-spin below. 



FIG. 4: The spatial distributions of the local density of states 
at the Fermi energy for the (0,0, z) surface in unit of c. For 
Na .75CoO 2 , (a) z = -0.0370, (b) z = 0.0, (c) z = 0.0462. 
For Na .2 5 CoO2, (d) z = -0.0370, (e) z = 0.0, (f) z = 0.0462. 



structures and the corresponding superconducting mech- 
anism, we strain the separation of the C0O2 layers for 
Nao.25Co02 and Nao.7sCo02 without taking into account 
the water molecules. The corresponding DOS with c = 
3.47A are displayed in Fig. 5. It is clearly seen that the 
hybridization between the cobalt and oxygen orbitals are 
considerably increased and the ferromagnetism prevails 
even in the high sodium content. Nao.25Co02, for exam- 
ple, the corresponding spin magnetization is increased 
to 0.785/x B , nearly two times of the corresponding value 
with normal separation of C0O2 layers. Thus, an effec- 
tive one-band model Hamiltonian may exist to describe 
the obtained electronic structure of this material. Actu- 
ally, the singularity in the DOS seen in Fig. 5 can induce 
strong quantum ferromagnetic fluctuations, from which 
an effective attraction among itinerant electrons is gener- 
ated as the same as the conventional low-temperature su- 
perconductors driven by the electron-phonon interaction 
|2.1 l24l I25I |2(| . The possible mechanism of superconduc- 
tivity in sodium cobalt oxides may thus originate from 
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FIG. 5: The density of states for Nao.7sCo02 and Nao.2sCo02 
with the enlarged separation of the C0O2 layers c = 3.47A. 
The density of states is in atomic unit (Ryd) -1 . The up-spin 
DOS is shown above the horizontal line and the down-spin 
below. 



In conclusion, pseudopotential method with local spin 
density functional approximation is used to calculate 
the electronic band structures of Na^OoC^ (x = 0.25, 
0.5, 0.75, and 1) with the structure relaxations. By 
increasing the Na content, both the hybridization be- 
tween cobalt and oxygen orbitals and the itinerant ferro- 
magnetism are decreased, and a phase transition is pre- 
dicted from a wide-band ferromagnetic to narrow band 
paramagnetic metals. It is the local distortions of the 
oxygen ions that reduces the itinerant ferromagnetism 
as increasing the sodium content. We also found that 
straining the C0O2 layers enhances both the hybridiza- 
tion and ferromagnetism. Whether the superconductiv- 
ity of Nao.35CoO2T.3H2O is related to a triplet (p-wave) 
electron pairing is required further investigation. 



the formation of a triplet (p-wave) electron pairing state. 
The similar mechanisms have been pro pos ed for the ma- 
terials UGe 2 , URhGe, and ZrZn 2 (Ref.[SllllI3)- 
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